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The structural and magnetic anomaly of the layered compound SrFe02 were examined by first 
principles density functional calculations and Monte Carlo simulations. The down-spin Fe 3d electron 
occupies the d^2 level rather than the degenerate {dxz, dy^) levels, which explains the absence of 
Jahn- Teller instability, the easy ab-plane magnetic anisotropy and the observed three-dimensional 
(0.5, 0.5, 0.5) antiferromagnetic order. Monte Carlo simulations show that the strong inter-layer 
spin exchange is essential for the high Neel temperature. 
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Perovskite oxides have attracted considerable inter- 
est due to their extensive applications in a number of 
technological areas. Among them, SrFeOa-a; and its re- 
lated iron perovskite oxides exhibit fast oxygen transport 
and high electron conductivity even at low temperatures 
[Hi 0, UI3, IB ■ It had previously been believed that the end 
member phases for these compounds are the orthorhom- 
bic brownmillerite SrFe02.5 {x = 0.5) and the cubic per- 
ovskite SrFeOs {x — 0). Very recently, the range of x 
was extended to a; = 1 by Tsujimoto et ai, who dis- 
covered that SrFe02 has planar Fe02 layers made up 
of corner-sharing Fe04 squares with high-spin Fe^+ (d^) 
ions, separated by Sr^+ ions d, 01 ■ 

SrFe02 exhibits interesting and apparently puzzling 
physical properties j^. First, if the lone down-spin elec- 
tron of a high-spin Fe^"*" {dP) ion at square-planar site 
occupies the degenerate {dxpdyz) orbitals, as expected 
by the crystal field theory [8|, SrFe02 should be sub- 
ject to orbital ordering or Jahn- Teller distortion when 
the temperature is lowered However, SrFe02 shows 
no structural instability and maintains the space group 
PA/mmm down to 4.2 K j^. Second, SrFe02 displays 
a three-dimensional (3D) antiferromagnetic (AFM) or- 
der with a very high Neel temperature (T/v — 473 K) 
0, which is even higher than that (^ 200 K) of FeO 
with a 3D structure. Such a high 3D AFM ordering 
temperature in a layered system is remarkable and un- 
expected, because T/v usually decreases drastically when 
the dimensionality decreases. Third, the powder neutron 
diffraction study shows that the magnetic moments are 
perpendicular to the c-axis (the local z-axis) 0, which 
is not consistent with the occupation of the (d^zjdyz) or- 
bitals with three electrons (loj . 

To probe the causes for these apparently anomalous 
structural and magnetic properties in SrFe02, we exam- 
ined the magnetic properties of SrFe02 by performing 
density functional theory (DFT) band structure and total 
energy calculations to evaluate its spin exchange interac- 



TABLE I: Spin exchange parameters (in meV) from the 
LDA-fU calculations. The spin exchange paths Ji, J2, J3, 
and J4 are defined in Fig. [T] Positive (negative) values indi- 
cate that the spin exchange interactions are AFM (FM). For 
SrFe02, both the experimental (expt.) and the optimized 
(opt.) crystal structures were used for the calculations. For 
CaFe02 and BaFe02, the optimized crystal structures were 
used. 



■h J2 J3 Ji 



SrFe02 (expt.) 7.04 2.18 0.43 -0.23 
SrFe02 (opt.) 7.91 2.29 0.30 -0.30 
CaFe02 (opt.) 8.90 3.24 0.35 -0.45 
BaFe02 (opt.) 5.81 1.34 -0.29 -0.16 




FIG. 1: (color online) Perspective view of the tetragonal 
structure of SrFe02. The large, middle, and small spheres 
represent the Sr, Fe, and O ions, respectively. The spin ex- 
change paths Ji, J2, J3, and J4 are also indicated. 



tions and performing Monte Carlo (MC) simulations to 
calculate the Neel temperature using the extracted spin 
exchange parameters. We show that the down-spin Fe 
3d electron in SrFe02 occupies the nondegenerate dz^ 
level rather than the degenerate {dxz, dyz) levels, which 
explains the absence of Jahn- Teller instability, the occur- 
rence of the easy a6-plane magnetic anisotropy as well as 
the strongly AFM coupling between the inter-layer and 
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FIG. 2: (color online) Electronic structure calculated for the 
FM state of SrFe02: (a) Dispersion relations of the up-spin 
and down-spin bands (solid and dashed lines, respectively) in 
the vicinity of the Fermi level, (b) Total DOS and PDOS 
plots for the Sr, Fe and O atom contributions, (c) PDOS 
plots for the Fe 3d orbitals, where the PDOS for the Fe d^o. 
orbital was highlighted by shading. 



intra-layer nearest-neighbor (NN) Fe^+ ions thereby giv- 
ing rise to the observed 3D (0.5, 0.5, 0.5) AFM order 
in SrFe02. MC simulations of specific heat show that 
the strong inter-layer spin exchange is essential for the 
high Neel temperature. Similar results are also found for 
the hypothetical isostructural compounds CaFe02 and 
BaFeOz. 

Our first-principles spin-polarized DFT calculations 
for MFe02 (M = Ca, Sr, Ba) were performed on the ba- 
sis of the projector augmented wave metho d ll lll encoded 
in the Vienna ab initio simulation package [12| using the 
local density approximation ^IS] and the plane-wave cut- 
off energy of 400 eV. To properly describe the strong 
electron correlation associated with the Fe Sd states, the 
LDA plus on-site repulsion U method (LDA-I-U) was em- 
ployed [3] . In the following, we report only those results 
obtained with V = 4.6 eV and J = eV on Fe [11], but 
the use of other U values between 3 — 6 leads to qualita- 
tively the same results. 

SrFe02 adopts the P4/mmm space group with a = 



3.985 A and c = 3.458 A at 10 



m spa 
K 



As shown in Fig. [1] 



the Fe02 layers separated by Sr^+ ions are stacked along 
the c axis. As a first step to discuss the magnetic proper- 
ties of SrFe02, the electronic structure of SrFe02 calcu- 
lated for its ferromagnetic (FM) state [q = (0.0, 0.0, 0.0)] 
is presented in Fig. [51 The band dispersion relations of 
Fig. mja) show that the FM state is metallic because the 
top portion of the up-spin valence bands (an O p and Fe 
d state) is above the bottom portion of the up-spin con- 
duction bands (a Fe p and Sr d state) at R (0.5, 0.5, 
0.5) point. The plots of the density of states (DOS) 
and the partial DOS (PDOS) in Fig. H^b) reveal that 
there is a strong hybridization between the O 2p and Fe 
3fi states in the valence bands. Further analysis shows 
that the largest coupling occur at R-point between Fe 
dr^2_y2 orbital and O Px,Py orbitals. When the PDOS 
plots of the Fe and O atoms in Fig. [D^b) are compared 
with those of the Fe 3d states in Fig. [2fc) , it is seen that 
the up-spin bands have Fe 3d energy levels below the O 
2p state, whereas the down-spin Fe 3d bands is above. 
This reflects the large exchange splitting of the Fe 3d 
levels, which is responsible for the high-spin state of the 
Fe^"*" ion. It is important to note from Fig. [5] that the 
occupied down-spin state has the character, not the 
doubly-degenerate {dxz, dyz) character as expected from 
crystal field theory for the Dj^h point symmetry 0, 0]. 
This is because in the layered structure, the energy of 
the c?22 state is significantly decreased due to the reduced 
Coloumb repulsion. Because the is non-degenerate, 
no Jahn- Teller type distortion is expected for SrFe02, as 
found experimentally Q. 

To determine the magnetic ground state and discuss 
the magnetic properties of SrFe02, we considered four 
more ordered spin states besides the FM state, namely, 
the AFl state with q = (0.5, 0.5, 0.5), the AF2 state with 
q = (0.0, 0.0, 0.5), the AF3 state with q = (0.5, 0.5, 0.0), 
and the AF4 state with q = (0.5,0.0,0.5). The experi- 
mentally observed AFM structure is AFl. Our LDA-I-U 
calculations show that all the AFM states are lower in en- 
ergy than the FM state, and the AFl state is the ground 
state, in good agreement with experiment [ij. These re- 
sults are consistent with the facts that in the high-spin 
(d^) configuration, there is no partially occupied d state 
to stabilize the FM phase The total DOS and PDOS 
plots presented in Fig. [3] show that SrFe02 in the AFl 
phase has a band gap as expected for this magnetic semi- 
conductor. Comparing to the electronic band structures 
of the FM phase, the band structure of the AFl phase 
exhibit an important difference; overall, the d bands are 
narrower in the AFl than in the FM state. In particu- 
lar, the up-spin d^^-y-^ band has a significantly narrower 
bandwidth in the AFl state. This observation is readily 
accounted for by considering a spin-1/2 square-net lattice 
with one magnetic orbital per site, the hopping integral 
t between adjacent sites and the on-site repulsion U . For 
the FM arrangement of the spins, the up-spin (or down- 
spin) states of adjacent sites are identical in energy and 
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FIG. 3: (color online) Electronic structure calculated for the 
AFl state of SrFe02: (a) Total DOS plot, which is identical 
for the up-spin and down-spin bands, (b) PDOS plots calcu- 
lated for the Fe 3d orbitals, where the shaded regions refer to 
the Fe 3d, 2 orbital contributions. 



the interaction energy between them is t, so the widths 
of the up-spin and down-spin bands are proportional to 
t. For the AFM arrangement of spins, the up-spin (or 
down- spin) states of adjacent sites differ in energy by U 
and the interaction energy between them is t'^ /U, so the 
widths of the up-spin and down-spin bands are propor- 
tional to t^/U. Since t ^ t'^/U for usual magnetic solids, 
the width of the electronic energy band is much wider for 
the FM state than for the AFM state. 

To extract the values of the four spin exchange param- 
eters Jl, J2 , J3 and J4 (see Fig. [1]) , we map the relative 
energies of the five ordered spin states (i.e., FM, AFl, 
AF2, AF3 and AF4) obtained from the LDA+U cal- 
culations onto the corresponding energies given by the 
Heisenberg spin Hamiltonian made up of the four spin 
exchange parameters. This mapping analysis was car- 
ried out as described in ref. 



17j . The obtained exchange 



parameters are summarized in Table. HI The intra-layer 
NN spin exchange Ji is quite strong, which is due to the 
strong 180° Fe-O-Fe superexchange between the dj;2_y2 
orbitals mediated by O Px and Py orbitals 18, 1^. The 
inter-layer NN spin exchange J2 is strongly AFM and is 
weaker than Ji only by a factor of ~ 3. The inter-layer 
spin exchange J2 originates from the direct through-space 
overlap between the dxz/dyz orbitals of Fe. 

To account for the observed anisotropic magnetic prop- 
erty of SrFe02, we carried out LDA+U calculations for 
the FM state by including spin-orbit coupling (SOC) in- 
teractions. These LDA-I-U-I-SOC calculations show that 
the state with the spin moments parallel to the a6-plane 
is more stable than the state with the spin moments par- 
allel to the c-axis by 4 meV/Fe, while the calculated or- 
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FIG. 4: (color online) Specific heat of SrFe02, C, calculated 
as a function of temperature T on the basis of the classical 
spin Hamiltonian defined in terms of the spin exchange and 
the spin anisotropy parameters. 



bital moment is 0.22 for the _L c spin arrangement, 
and 0.01 /is for the || c spin arrangement. This easy 
a6-plane anisotropy, which is in accord with experiment 
g|, can be explained by analyzing the SOC Hamiltonian 
20|. Since the up-spin and down-spin d bands are well 
separated due to the large exchange splitting, one can 
neglect interactions between the up-spin and down-spin 
states under the SOC. With 9 and as the zenith and 
azimuth angles of the magnetization in the direction n{9, 
(f)), the spin-conserving term of the AL-S operator (A < 
for Fe2+) is given by [2l| 



XSn{L^ cose+ -L+e-"^ sin( 



1 



L_e"^sin6i) (1) 



If SrFe02 were to have one down-spin electron in the 
doubly-degenerate level [dxz,dyz), it should have uniax- 
ial magnetic properties with the spin moments parallel 
to the c-axis {9 = 0°) according to the degenerate per- 
turbation theory [l^]. Our calculations show that, for 
the down-spin part, the highest occupied Fe 3d level is 
dz2 while the lowest unoccupied Fe 3d level is {dxz, dyz). 
When the spin lies in the a6-plane (9 = 90°), the mixing 
between d^^ and dxzjdyz due to the raising and lowering 
operators is the largest. Thus, nondegenerate pertur- 
bation theory shows that SrFe02 has an easy ab- plane 
anisotropy with a relative large orbital moment. 

To explain why the layered compound SrFe02 has a 
very high Neel temperature (473 K), we perform MC 
simulations for a 12 x 12 x 12 supercell based on the 
classical spin Hamiltonian: 



2 

iz ' 



(2) 



<^J> 



where the spin exchange parameters Jij are those defined 
in Fig. m £1 = 1 meV (i.e., DSf, ^ AD = A meV) is the 
spin anisotropy parameter, and S — 2. To obtain T/v, 
we first calculate the specific heat C = {{E^) - {E)^)/T^ 
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after the system reaches equihbrium at a given tempera- 
ture (T) . Then T/y can be obtained by locating the peak 
position in the C{T) vs. T plot, shown in Fig. 21 For 
SrFe02 the calculated Tjv is 354 K, which is in reason- 
able agreement with the experimental value of 473 K. 
We should note that a smaller U value leads to larger 
exchange parameters Ji since Ji oit^ /U , and J2 changes 
a little, thus Tjv will be higher. For instance, \iU = 4.0, 
then Ji — 8.68 meV and J2 — 2.23 meV, resulting in 
Tn = 453 K. 

It was shown that large magnetic anisotropy en- 
ergy can stabilize long-range magnetic order in a one- 
dimensional (ID) monatomic Co chains .22]. To deter- 
mine how Tjv is affected by the different spin exchange 
parameters and the spin anisotropy, we performed three 
additional simulations; one with the spin anisotropy D 
neglected, one with the inter-layer exchange parameters 
(J2 and J4) neglected, and one with both the inter- 
layer spin exchange parameters ( J2 and J4) and the spin 
anisotropy D neglected. The resulting C(T) vs. T plots 
are presented in Fig. (4] As can be seen, when D is dis- 
regarded, T/v is only slightly lower (347 K). If only the 
inter-layer exchange parameters are neglected, the C{T) 
vs. T plot shows a broad peak at 220 K. A broad peak oc- 
curs at 208 K if both the inter-layer exchange parameters 
and the spin anisotropy are neglected. At any nonzero 
temperature, a ID or 2D isotropic Heisenberg model with 
finite-range exchange interaction can be neither FM nor 
AFM 123|. Thus, the broad peak in the C{T) vs. T plot, 
obtained when the inter-layer spin exchange interactions 
are neglected, indicates the presence of short-range or- 
der. Our MC simulations indicate that the inter-layer 
interactions are primarily responsible for the high Tjv of 
SrFe02. 

It should be noted that CaFeOa-x and BaFeOa-x have 
structures and properties similar to those of SrFeOa-x- 
Thus, it is of interest to consider the structural and mag- 
netic properties of hypothetical CaFe02 and BaFe02 as- 
suming that they are isostructural with SrFe02. For this 
purpose, the structures of MFe02 (M = Sr, Ca, Ba) 
were optimized by performing LDA-I-U calculations for 
the AFl state. The optimized lattice constants a and c 
for SrFe02 are 3.92 A and 3.40 A, respectively, which are 
in good agreement with experiment. The calculated lat- 
tice constants are a = 3.86 A and c — 3.12 A for CaFe02, 
and a = 3.98 A and c = 3.79 A for BaFe02. The trend in 
the lattice constants is consistent with the ionic radii of 
Ca^+, Sr^+, and Ba^+. The exchange parameters calcu- 
lated for the optimized MFe02 (M — Sr, Ca, Ba) struc- 
tures are listed in Table. HI All these compounds should 
have the AFl state as the ground state since the spin 
exchange interactions are dominated by AFM Ji and J2 . 
Among MFe02 (M = Ca, Sr, Ba), CaFe02 has the largest 
Ji and J2 values, while BaFe02 has the smallest Ji and 



J2 values. This trend reflects the fact that the strengths 
of these interactions increase with decreasing the lattice 
constants. Thus, CaFe02 is predicted to have a higher 
Tjv than does SrFe02. 

In summary, SrFe02 has no Jahn- Teller instability be- 
cause the occupied down-spin d-level is (1^2 . The Neel 
temperature of SrFe02 is high because the intra-layer NN 
spin exchange is strong while the inter-layer NN spin ex- 
change is substantial. The in-plane magnetic anisotropy 
of SrFe02 arises from the SOC-induced interaction be- 
tween the (1^2 and {dxz, dyz) states. 
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